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^ I Abstract. I explain the Sato-Lee (SL) model and its extension to the neutrino-induced pion produc- 

tion off the nucleon. Then I discuss applications of the SL model to incoherent and coherent pion 
CN ' productions in the neutrino-nucleus scattering. I mention a further extension of this approach with 

a dynamical coupled-channels model developed in Excited Baryon Analysis Center of JLab. 
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3 : INTRODUCTION 

Neutrino oscillation experiments have been actively conducted in the last decade, and 
^ , will be so in the forthcoming decade. Because those experiments detect the neutrino 

■^ I through the neutrino-nucleus (v-A) scattering, understanding of the v-A scattering is a 

^ ' prerequisite for a successful interpretation of data. Some of the experiments measure the 

^ . neutrino in the energy region of sub- and few-GeV where dominant processes are quasi- 

elastic nucleon knockout (QE) and single pion (In) production via the A-excitation. For 
O ' QE, the elementary amplitude is reasonably well-known, and the challenge is to incor- 

porate the nuclear correlation in the initial state and the final state interaction. Although 
there has been a reasonable success in describing QE in the electron-nucleus scattering, 
it was reported that the same framework does not work well for the v-A scattering [1]. 
^ ■ From here, I focus on the In productions in v-A scattering in the A region, which con- 

stitute the dominant background in the neutrino oscillation experiments. In addition to 
the difficult problem of the nuclear effects, relevant elementary amplitudes, or dynam- 
ical models which generate them, also have to be carefully studied. Those dynamical 
models are developed through a careful analysis of data for electroweak 1 n production 
off the nucleon. Actually, there have been active developments of such dynamical mod- 
els, motivated by extensive experiments of photo- and electro meson-productions in the 
resonance region. These experiments aim to test resonance properties predicted by QCD- 
inspired models and Lattice QCD. A dynamical model for Ik production developed in 
this way provides a good starting point to study the neutrino-induced In production off 
the nucleon, because of the close relation between the weak and the electromagnetic 
currents. Furthermore, the dynamical model offers a good basis to study In production 
in the v-A scattering. 

Thus first, I give a brief description of a dynamical model, i.e., the Sato-Lee (SL) 



model [2], for the In photo-production off the nucleon. (For a fuller discussion, consult 
Ref. [2].) Then I discuss the extension of the SL model to the weak sector, done in 
Ref. [3]. The elementary amplitudes generated by the SL model has been applied to 
pion productions in v-A scattering. I describe the work done in Ref. [4] where the authors 
studied the quasi-free A-excitation followed by the single pion production. I also discuss 
the coherent pion production off a nucleus studied with the SL model [5]. Finally, I 
discuss a possible future development. 

SATO-LEE (SL) MODEL 

In the SL model, one starts with a set of phenomenological Lagrangians, and derive an 
effective Hamiltonian using a unitary transformation. The effective Hamiltonian for pion 
photoproduction can be written as follows: 

^eff = Ho + VytM + VyK+Tj^pj^A + ^jN^Ai (1) 

where Hq is the free Hamiltonian, Vjzf^ and Vyy^ are respectively non-resonant nN — )■ ;rA'^ 
and yN -^ kN potentials, and are composed by the Born diagrams, ^-channel p and (O 
exchange terms, and the crossed A term. Bare vertices for kN f-> A and yA^ ^ A tran- 
sitions are respectively denoted by T^zn^a and TyM<r^A- With the effective Hamiltonian, 
we can derive unitary pion photoproduction amplitude as 



^h^KN{E)TyM^A{E) 

E-niA-T^AiE) 



TyAE) = tyAE) + ^r^r^ ; :^^ (2) 



where the first (second) term is the nonresonant (resonant) amplitude, and E is the total 
energy of the pion and nucleon; m./^ is the A bare mass. The nonresonant amplitude is 
calculated by 

tyn{E) = Vyjt + tjtN{E)GnN{E)vyn, (3) 

where G;;:yv is the nN free propagator, and tjj:pj is obtained by solving Lippmann- 
Schwinger equation which includes Vjtf^. In Eq. (2), the A vertices are dressed as 

^yN^A{E) = Ty]^^i^ + VyTiGjiN{E)T^i^_^i^{E), (4) 

^a-^kn{E) = [l+tji]\i{E)G,tN{E)]TA^j^f^ , (5) 

so that dynamical pion cloud effect is taken into account as a consequence of the 
unitarity. The A self-energy in Eq. (2) is given by 

'^a{E) =T^j^^/i^G^N{E)r/^_^j^j^{E). (6) 

The ttN scattering amplitude is calculated similarly. Thus one first determine the strong 
interactions by analyzing nN data in the A-region. Then adjustable parameters rele- 
vant to electromagnetic interactions are fixed by analyzing yN -^ kN data. With this 
approach, the SL model has been shown to give a reasonable description for the kN, 
yN — )> TtN [2] and eN -^ e'nN reactions [6] in the A-region. 
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FIGURE 1. Total cross sections for v^A^ -^ pL nN. For an explanation of each curve, see the text. Data 
are fromRef. [7]. 



iTT PRODUCTION IN NEUTRINO-NUCLEON SCATTERING 

With the SL model for pion photo-(and electro-) production discussed above, it is 
straightforward to extend it to the weak sector, as has been done in Ref. [3]. One just 
need to replace the electromagnetic current with the with V^ —^^ where V^ (A^) is the 
weak vector (axial- vector) current. The vector current conservation hypothesis tells us 
that the weak vector current is obtained from the isovector part of the electromagnetic 
current by the isospin rotation. Thus the remaining part is the axial-current. In Ref. [3], 
the authors parametrized the least known axial-vector A^A transition matrix element, 
more specifically form factors, as 



gAN^{Q^)=gANmR{Q^)GA{Q^ 



(7) 



where Ga{Q^) = \/{\+Q^/Mlf with Ma = 1.02 GeV. The coupling ^awa(O) is related 
to the nucleon axial coupling gA using the nonrelativistic constituent quark model. The 

remaining correction factor, R{Q^) = (1 -\- aQ^)e^^^ , is assumed to be the same as 
that used for the yN -^ A form factor, and thus is determined by analyzing the pion 
electroproduction data. 

The total cross sections predicted by this model is compared with data in Fig. 1. It 
turns out that the full calculation (solid curves) shows a good consistency with the data. 
If we turn off the meson cloud effect, then we obtain the dotted curves, indicating the 
significant effect. We further turn off the contribution from the bare A^A transition, then 
we obtain the dashed curve, showing the non-resonant contribution. Although the non- 
resonant contribution is smaller than the resonant one, it is still important to get a good 
agreement with data because it can interfere with the resonant amplitude. 



In PRODUCTION IN NEUTRINO-NUCLEUS SCATTERING 



The SL model discussed in the previous section can be applied to the neutrino-nucleus 
interaction in the A-region, which has been conducted in Ref. [4] for the ^^C target. A 
unique feature of this application is that the SL model treats resonant and non-resonant 



> 

C5 



o 

I 



0.8 


'"s 




, ' 6=10" 


0.6 


/ '. 




/ /-"H 


0.4 


■ /'//■^■■ 




/■///^ 


0.2 
n 





> 

C5 

c 

I 



H-UU 

300 


■ 


A 


200 




\fe^^^%«' 




57 


Vy/ \ 


100 
n 


! / 


/\ 



1.1 



1.2 1.3 

W [GeV] 



1.4 



0.8 



1 1.2 

W [GeV] 



1.4 



FIGURE 2. (Left) Nuclear effects on differential cross sections for Ve +'^C — > e^%X at = 10°, 
normalized with the target mass number (Right) Differential cross sections for e^ +'^C — ^ e^X at 
Eg = 1.1 GeV with Qf, = 37.5°. Data are from Ref. [9]. For a description of each curve, see the text. 



mechanism on the same footing so that the amplitude is unitary, while most previous 
works considered only resonant mechanisms. The challenge here is to incorporate el- 
ementary amplitudes generated by the SL model with various nuclear effects such as: 
the nuclear correlation effect in the initial state; the Pauli blocking on the final nucle- 
ons; the final state interactions (including pion absorption); the medium effect on the 
A-propagation. The authors of Ref. [4] considered the initial nuclear correlation using 
the spectral function [8], and the Pauli blocking using the Fermi gas model. 

Left panel of Fig. 2 shows the nuclear effect on the differential cross sections for 
Vg+^^C —> e^nX, normalized with the mass number, at the lepton scattering angle 
6 = 10° as a function of the final ^rA'^ invariant mass W. The dashed curve shows the 
differential cross sections for the neutrino-induced pion production off the free nucleon, 
averaged over the free proton and neutron. With the Fermi gas effect on the initial 
nucleon distribution, we obtain the dashed-double-dotted curve. We can see that the 
Fermi motion broaden the A-peak. By considering the Pauli blocking in addition to the 
Fermi motion, we obtain the dash-dotted curve. The Pauli blocking reduces the forward 
cross sections by about 20%. Finally, the solid curve is obtained by replacing the Fermi 
gas model with the spectral function taken from Ref. [8]. The spectral function further 
broadens the peak, and reduces the height of it by about 20%. 

For demonstrating the validity of the approach, the model prediction for e^ +^^C — ^ 
e X is compared with data in Fig. 2 (right). The QE and A peaks reasonably reproduce 
data. However, the dip-region is rather underestimated, indicating the need of going 
beyond the impulse approximation, and/or more elaborate treatment of the nuclear 
correlation. 



COHERENT k PRODUCTION 

The SL model has been applied to the coherent pion production on ^^C in Ref. [5]. The 
approach taken in Ref. [5] is to combine the elementary amplitudes from the SL model 
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FIGURE 3. (Left) Differential cross sections for the coherent pion production off ^^C. Data are from 
Ref. [10]. (Right) The pion momentum distribution for the charged-current coherent pion production in 



V— C scattering. For a description of each curve, see the text. 



with the A-hole model. This approach allows us not only to implement the nuclear effects 
such as modification of the A-propagation and the pion absorption, but also to describe 
n-A scattering, coherent tt photoproduction and coherent n production in v-A scattering 
in a unified manner. Thus, we can fix parameters relevant to the medium modification 
of the A-propagation by analyzing n-A (total and elastic) scattering, and then we can 
predict the coherent pion productions. 

Figure 3 shows the nuclear effects and the predictive power of the model. The dashed 
curves in the left and right panels include neither the medium modification on the A- 
propagation nor the final state interaction between the pion and nucleus. The shape of 
the curves are determined by the elementary amplitude, the nuclear form factor, and 
the phase-space factor. With the medium effect on the A, we obtain the dotted curves. 
By further turning on the final state interaction, we obtain the solid curves. The nuclear 
effects are very significant, and bring the calculation into a good agreement with data 
for the photo-production. This is an important test of the model. 

The dash-dotted curves are obtained by turning off the non-resonant amplitudes. 
By observing the difference between the solid and dash-dotted curves, we can see a 
significant contribution from the non-resonant amplitude, even in the A-region. This is 
in contrast with the finding in Ref. [11] that the non-resonant amplitude plays essentially 
no role. It is noted that Ref. [11] used a tree-level elementary amplitude, while Ref. [5] 
used a unitary one from the SL model. The difference in the reaction mechanism may 
be responsible for differences observed in theoretical predictions of the pion momentum 
distributions and Ey -dependence of the total cross sections. 

We can average the total cross sections using the neutrino flux from experiments. For 
the charged-current (CC) process, we use the flux from K2K, and obtain 6.3 x lO^^'^cm^ 
which is consistent with the report from K2K [12], < 7.7 x lO^^'^cm^. For the neutral- 
current (NC) process, we use the flux from MiniBooNE to obtain 2.8 x lO^^^cm^ 
which is still consistent within the rather large error bar of the preliminary report [13]: 
7.7 ± 1.6 ± 3.6 X lO^^^cm^. However, the CC/NC ratio is not in agreement with the 
recent report [14], as no theoretical calculations are not. 



FUTURE DEVELOPMENT 

Having seen reasonable descriptions of pion productions in neutrino(photon, electron)- 
nucleus scattering with the SL model plus nuclear effects, it is highly hoped to extend 
this approach to higher mass resonance region. This is because a model that covers 
the region from the A to DIS is very useful for neutrino experiments. In this energy 
region, several hadronic channels couple, and 2k production reactions occupy quite a 
little portion of final states. Thus we need a dynamical model that takes care of channel- 
couplings, and treat the single and double meson productions on the same footing. In 
this context, continuous effort made at the Excited Baryon Analysis Center (EBAC) in 
JLab [15] is quite encouraging. The EBAC has been analyzing world data of yN, kN — )■ 
kN^ tzkN, r\N, KY reactions in the resonance region with a dynamical coupled-channels 
model (EBAC-DCC model), and aim to extract resonance information. The EBAC-DCC 
model is an extension of the SL model by extending the coupled-channels from kN 
to KN^r]N^TZKN{Tt/S.^oN^pN),KA,K'L, and also by including higher resonance states. 
It has been demonstrated that the EBAC-DCC model gives a reasonable description 
of pion- and photo-induced meson production reactions from the A to higher mass 
resonance region [16]. An extension of the EBAC-DCC model to the weak sector 
and neutrino-nucleus reaction, as done with the SL model, seems a promising future 
direction. 
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